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Abstract 

We present a relation which connects the propagator in the radial 
(Fock-Schwinger) gauge with a gauge invariant Wilson loop. It is closely 
related to the well-known field strength formula and can be used to cal- 
culate the radial gauge propagator. The result is shown to diverge in 
four-dimensional space even for free fields, its singular nature is however 
naturally explained using the renormalization properties of Wilson loops 
with cusps and self-intersections. Using this observation we provide a 
consistent regularization scheme to facilitate loop calculations. Finally 
we compare our results with previous approaches to derive a propagator 
in Fock-Schwinger gauge. 
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1 Introduction 



While perturbation theory for gauge fields formulated in covariant gauges is 
very well established ||l| many aspects of non-covariant gauges are still under 
discussion. In principle one expects physical quantities to be independent of 
the chosen gauge. However this might lead to the naive conclusion that a quan- 
tum theory in an arbitrary gauge is simply obtained by inserting the respective 
gauge fixing term and the appropriate Faddeev-Popov ghosts in the path inte- 
gral representation and reading off the Feynman rules. Unfortunately it is not 
so easy to obtain the correct Feynman rules, i.e. a set of rules yielding the same 
results for observable quantities as calculations in covariant gauges. Prominent 
examples are formulations in temporal and axial gauges. Such gauge choices 
are considered since one expects the Faddeev-Popov ghosts to decouple. How- 
ever problems even start with the determination of the appropriate free gauge 
propagators. Temporal and axial gauge choices yield propagators plagued by 
gauge poles in their momentum space representation. These are caused by the 
fact that such gauge conditions are insufficient to completely remove the gauge 
degrees of freedom. The correct treatment of such poles can cause ghost fields 
to reappear can break translational invariance |^ or both Q. While these 
problems seem to be "restricted" to the evaluation of the correct gauge propaga- 
tors and ghost fields, the necessity of introducing even new multi-gluon vertices 
appears in the Coulomb gauge |^ . These additional vertices are due to operator 
ordering problems which are difficult to handle in the familiar path integral ap- 
proach. They give rise to anomalous interaction terms at the two-loop level 
and cause still unsolved problems with renormalization at the three-loop level 

0- 

In this article we are interested in the radial (Fock-Schwinger) gauge condition 

Xf,A^'{x) = 0. (1.1) 

It found widespread use in the context of QCD sum-rules (e.g. There it is 
used as being more or less synonymous to the important field strength formula 



^-^(a;) = Jdssx-'F^^isx) (1.2) 



which enormously simplifies the task of organizing the operator product ex- 
pansion of QCD n-point functions in terms of gauge invariant quantities by 
expressing the gauge potential via the gauge covariant field strength tensor. It 
was introduced long ago [||, ||l^ and rediscovered several times (e.g. |ll|]). 

Only a few efforts have been made to establish perturbation theory for radial 
gauge. The main reason for this is that the gauge condition breaks translational 



invariance since the origin (in general an arbitrary but fixed point z, c.f. (1.6)) 
is singled out by the gauge condition. Thus perturbation theory cannot be 
formulated in momentum space as usual but must be set up in coordinate space. 
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1 Introduction 



The first attempt to evaluate the free radial propagator was performed in Q . 
Later however the function T^i,{x, y) presented there was shown to be not sym- 
metric Moreover it could not be symmetrized by adding Tv^{y,x) since 
the latter is not a solution of the free Dyson equation. It was even suspected 
in ||l3| that it might be impossible to find a symmetric solution of this equation 
in four-dimensional space, due to the appearance of divergences even on the 
level of the free propagator when one uses the field strength formula to derive 
a free propagator. Indeed we agree with this statement in principle, but we 
will present an explanation for this problem and a way to bypass it. Other 
approaches to define a radial gauge propagator try to circumvent the problem 
(e.g. by sacrificing the field strength formula as given in ( |l.2| ) which was 
one of the main reasons the gauge became popular in non-perturbative QCD 
sum rule calculations in the first place. If we are not prepared to do so we 
are forced to understand the origin of the divergences that plague most of the 
attempts to define even free propagators in radial gauges and see whether they 
can be dealt with in a satisfying manner. 

In Section we will make the first and decisive step in this direction by explor- 



ing the completeness of the gauge condition (1.1) and its relation to the field 
strength formula and developing a new representation of the gauge potentials 
via link operators. 

In Section ^ we use this information to relate the divergences encountered in 
some of the attempts to define radial propagators to the renormalization prop- 
erties of link operators. We find that even free propagators in radial gauge 
may feel remnants of the renormalization properties of closed, gauge invariant 
Wilson loops. Surprising as this seems to be superficially it is not impossible 
however if we recall that the inhomogeneous term in the gauge transformation 
has an explicit 1/g factor in it. As a result we are able to define a regularized 
radial propagator using the field strength formula and established regularization 
procedures for link operators. 

Section ^ will be devoted to demonstrate the consistency of our approach by 
calculating a closed Wilson loop using our propagator and relating the steps to 
the equivalent calculation in Feynman gauge. 

In Section ^ we obtain an explicitly finite version of our propagator by complet- 
ing the renormalization program developed for link operators before we sum- 
marize and compare our results to other approaches in the literature in ^ and 
shortly discuss the next steps in the program of establishing a new perturbative 
framework in radial gauges which - although the steps to be performed are quite 
straightforward - we will postpone for a future publication. 

In the following we work in a ZJ-dimensional Euclidean space. The vector po- 
tentials are given by 

A,{x)^Al{x)t, (1.3) 
where ta denotes the generators of an SU{N) group in the fundamental repre- 
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sentation obeying 

[ta,tb]^ifabct^ (1.4) 

and 

t'citatb) = ^Sab- (1.5) 

In general the radial gauge condition with respect to z reads 

(x- z)^A^ix) = 0. (1.6) 

For simplicity we take z = after Section ^. The results nevertheless can be 
easily generalized to arbitrary values of z. 



2 The gauge condition revisited 



Before we can go ahead and tackle the problem of divergences in the radial 
gauge propagator we have to establish a clearer picture of the uniqueness of the 
gauge condition we are about to implement. After all, if we do not succeed to 
fix the gauge completely we might be naturally confronted with divergences - 
if not at the free level then later in perturbative calculations. They would be 
a simple consequence of the incompleteness of the gauge fixing and the zero 
modes of the propagator which would then necessarily be present. This point 
has caused a still continuing discussion for the case of axial gauges (e.g. Q) but 
is only briefly mentioned in the context of radial gauges (e.g. p^). 

Readers who are not interested in the discussion of (in)completeness of radial 
gauge conditions might skip the following considerations without getting lost 
and start reading again after eq. (2.6). 



To clarify the question whether the gauge condition (1.6) is sufficient to com 



pletely fix the gauge degrees of freedom we have to catalogue the gauge trans- 
formations U\B]{x ) w hich transform an arbitrary vector potential B into the 
field A satisfying (L£). A gauge condition is complete if U[B]{x) is uniquely 
determined up to a global gauge transformation. In other words, we want to 
find all solutions of 



(x-z)^ U[B]{x) 



Bi'ix) 



1 



U[B]-^{x) = 0. 



(2.1) 



It is easily checked that we have an infinite family of such solutions which can 
all be cast in the form of a product of two gauge transformations of the form 



Here 



U[B]{x) = V{z{x))U[B]{z{x),x) . 



U[B]{z{x),x) ^Vexpig I duj^,B^'{uj) 



(2.2) 
(2.3) 
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2 The gauge condition revisited 





a) 



b) 



Fig. 1: a) straight line path in tire linlcs for two points x and y. 
spherical hyper-surface z{x). 



b) Example for a 



is a link operator whose geometric ingredients are parameterized via its end- 
points X and z(x) and the straight hne path co between them, V denotes path 
ordering. 

In particular z{x) is the point where a straight line from z through x and a given 
closed hyper-surface around z intersect. Since there is a unique relation between 
these points and the hyper-surface we will also refer to the hyper-surface itself 
by z{x). This geometry is illustrated in Fig. |l|. Both the detailed form of z{x) 
and the local gauge transformation V{x) are completely unconstrained as long 
as {x — z).d^z{x) = 0. In short, they parameterize the residual gauge free- 
dom not e limi nated by {lA). Note that while V{x) is completely arbitrary the 
solutions ( ^.2[ ) ask only for its behavior at the given hyper-surface z{x). The 
simplest and most intuitive choice for z(x) is a spherical hyper-surface around 
z. Introducing the appropriate spherical coordinates it becomes obvious that 
V{z{x)) parameterizes gauge transformations which purely depend on the an- 
gles. Clearly the radial gauge co ndit ion (1.6) cannot fix the angular dependence 
of any gauge transformation in (2.1). 



To eliminate t he r esidual gauge freedom one has to impose a condition which is 
stronger than (1.6) and suffices to pin down V{x) up to a global transformation. 
A possible choice for such a gauge fixing would be the condition 




1 



du;.Aiu;)+ dSj -{z{x),y) 9^.A(y) =0 



(2.4) 
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which in addition to the vanishing of the radial component of the gauge potential 
also implements a covariant gauge on the hyper-surface z{x). Such a gauge for 
arbitrary z{x) would immediately force us to introduce ghosts into the path 
integral. Moreover the field strength formula would also be lost as we will 
illustrate below. 

There is one exception to these unwanted modifications however, which may 
be implemented by contracting the closed surface z(x) to the point z. Then 
the influence of V{x) becomes degenerate with a global transformation and the 
gauge is completely fixed. Incidentally this is also the only case which entails 
the field strength formula. To see this we use 

SU{x, z) = ig \ Af,{x)U{x, z)dx^' - U{x, z)A^{z)dz''+ 



ds [U{x,Wx)Fi_^^{wj;)U{wj;, z)] 



dwit ( dwZ 



ds V rfa;" 



-dx" 



^dz- 
dz°' 



(2.5) 

(see e.g. [0) to differentiate the link operators in the gauge transformation 
(^.2[) in order to find an expression for the radial gauge field: 



A^;%x) = U[A]iz,x) 
1 



1-9 



U[A]{x,z) 



dio^ 

ds s—U[A]{z,oj) F,^{uj)U[A]{uj, z) 



(2.6) 



This is nothing but (1.2) for arbitrary z (note that in this case oj — lo{s) is 
simply given hy lo{s) — z + (x — z)s.) This simple result is only true since 
dfj,z{x) = dfj,z = 0. For general z{x) there would be an additional term in the 
above formula reflecting the residual gauge freedom encoded in V{z{x)). 

This sets the stage for a further exploration of the radial gauge in a context 
where we can be sure of having completely flxed the gauge in such a way that 
the field strength formula is guaranteed to be valid. Before we go on to study- 
ing the consequences the above has for the implementation of propagators we 
will introduce yet another representation of the gauge field in this particular 
complete radial gauge - this time solely in terms of link operators. 

From now on we will assume the reference point z to be the origin, but it will 
always be straightforward to recover the general case without any ambiguities. 
We will also suppress the explicit functional dependence of link operators on 
the gauge potential A for brevity. 
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3 The Radial Gauge Propagator 



Let us start with a link operator along a straight line path 



U{x,x') = Pexp 



(2.7) 



where now w(s) := x' + s{x — x') . According to (2.5) we have 



d-^U{x,x')=ig 



A^{x) - I ds s —J— U{x,uj) F^^{lu)U{uj,x) 



U{x,x') (2.i 



which can be used to express the vector potential in terms of the link operator 

lim U{x,x') = igA^(x) . (2.9) 

x' — >x 

In the case at hand the fact that U{0, a;) = 1 in any of the x.A{x) = gauges 
allows us to introduce a new gauge covariant representation 



A'^'^ix) = - lim df, [U(0, x) U(x, x') U(x', 0)1 



(2.10) 



for the Fock-Schwinger gauge field. It is easy to see that this is indeed equivalent 
to the field strength formula as given in_(1.2) and consequently satisfies the same 
complete gauge fixing condition (i.e. (2^) for z{x) z): 



AJf^(x) = - lim dl [[/(O, x) U{x, x') U{x', 0)] 
ig x'^x ^ 

= — lim \d''U{0,x)U{x,x')U{x',0) 
ig x'^x ^ 



'U{Q,x)dlU{x,x')U{x\Q)\ 



1 



= —dl U(0, x) U(x, 0) + U{0, x) AJx) U{x, 0) 
1 



ds sx^'Flfisx) 



(2.11) 



where the last step uses (p.§|), mirroring the relations in (2.6) for z = 



3 The Radial Gauge Propagator 

Having established the complete gauge fixing we are interested in, it is now 
straightforward to devise expressions for the propagator as a two-point function. 
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According the above we know that 

= hm d-^dy (C/(0, x) U{x, x') U{x', 0) ® f/(0, y) U{y, y') U{y', 0)) 

y' —*y 
1 1 

ds Jdtsx'^ty^ (U{0,sx)F^f,{sx)U{sx,0) (^U{0,ty)Ff3^{ty)U{ty,0)) . 



(3.1) 

Since we are in a fixed gauge it makes sense to perform a muhiplet decomposition 
and for instance extract the singlet part of this propagator. The latter reduces 
to the free propagator in the limit g —^ 0. 

We define 

A[2 _ 1 

tr {A^{x)A,{y)) = tr(t,4) {Al{x)At{y)f^^"^' = -—- D^,{x, y) (3.2) 



= :5-'''D^,{x,y) 

to extract 

{Al{x)Al{y))f:t" = ^'''j^^' {A,{x)A,{y)Ua - '^"'a^ ^^'^^ 
X -i- hm d^dl tr (C/(0, x) U(x, x') U{x' , 0) C/(0, y) U{y, y') U{y', 0)} . 

V —'V 

Obviously 

W^{x, x', y') := 1 tr ([/(O, x) U{x, x') U{x', 0) U{0, y) U{y, y') U{y', 0)) (3.4) 
is a gauge invariant Wilson loop. Its geometry is illustrated in Fig. ||. 



On the other hand, using the second expression in ( |3.l[ ) we have an equiva- 
lent representation for the singlet part of radial gauge propagator via the field 
strength formula: 

(^;:w4(2/))r/''* = '5'^'^^ (3.5) 

1 1 

X Jds Jdtsx" tyf^ ti {U{0, sx) Faf,{sx) U{sx, 0) U{0, ty) F0^{ty) U{ty, 0)) . 



Modanese tried to calculate the free radial gauge propagator from ( |3.5D in 
a D dimensional space-time.^ Unfortunately one gets a result which diverges in 
the limit D A. 



'In fact he discussed the Abelian case but this makes no difference for free fields. 
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3 The Radial Gauge Propagator 




Fig. 2: The Wilson loop Wi which is int imately connected with the 
radial gauge propagator according to ({3.i), (3.i.)) 



Since we have performed a complete gauge fixing (at least on the classical level), 
this comes as a surprise since we certainly do not expect zero mode problems 
to come into the game as a possible explanation and consequently a way out. 
Does this mean we are trapped at a dead end or is there another explanation 
for this seemingly devastating discovery? 

Before we try to answer this question let us briefly recapitulate how this di- 
vergence makes its appearance: Since the right hand side of ( |3.5| ) is gauge 
invariant we can choose an arbitrary gauge to calculate it. For simplicity we 
take the Feynman gauge with its free propagator 



T{DI2 - 1) 
47r^/2 



_„^2^1-D/2 



(3.6) 

Using the free field relations U{a, b) — \ and F^^i, = df^Ai, — d^A^^ we get (for 
more details see Appendix 0) 



{Al{x)Al{y))l 



rad 

-^-^i!^^'ids}dtsx-ty^ 



1 1 

-dl fdsx.iisx-yfj'-^/'-dy fdty^lix-tyf]'-^/' 



2ll-_D/2 
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1 1 

+ d^^dy Jds Jdtx-y[{sx-tyf]'-°/^ ]. (3.7) 





1 

Thus the radial gauge propagator is singular for arbitrary arguments x and y 
— with one remarkable exception: It is easy to see that it vanishes for a; = or 
y = 0. This is simply a consequence of our task to preserve the field strength 



formula (1.2) which forces the vector field to vanish at the origin (in general at 



the reference point z). 

The observation that the radial gauge propagator as calculated here diverges 
in four-dimensional space raises the question, whether it is perhaps impossible 
to formulate a quantum theory in radial gauge. This would suggest that the 
radial gauge condition - in the form that facilitates the field strength formula 
- is inherently inconsistent ( "unphysical" ) in contrast to the general belief that 
it is "very physical" since it allows to express gauge variant quantities like the 
vector potential in terms of gauge invariant ones. To answer this question we 
have to understand where this divergence comes from. In the following we will 
see t hat for this purpose the complicated looking Wilso n lo op representation 
( ^.3[ ) is much more useful than the field strength formula ( |3.5|) . (Note, however. 



that the result for the free propagator (3.7) of course will be the same 



It is well-known that Wilson loops need renormalization to make them well- 
defined (see e.g. and references therein). The expansion of an arbitrary 
Wilson loop 

W{C) = ^tr (vexp 



N 

in powers of the coupling constant is given by 



ig f dx'^Af^{x) 



c 



(3.8) 



x<dc{xi > ■ ■■ > Xn) trG^j...^„(a;i, . . . ,x„) (3.9) 
where Qci^i > ■ ■ ■ > Xn) orders the points Xi, . . . , x„ along the contour C and 
Gf^^...f^^{xi, . . . ,x„) {Af^^{xi) ■ • •yl^„(x„)) (3.10) 
are the Green functions. 



In general Wilson loops show ultraviolet singularities in any order of the coupling 
constant. If the contour C is smooth (i.e. differentiable) and simple (i.e. without 
self-intersections) the conventional charge and wave-function renormalization — 
denoted by TZ in the following — is sufficient to make W{C) finite. We refer 
to p9[ | for more details about renormalization of regular (smooth and simple) 
loops. 
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In our example we must apply the renormalization operation TZ to Wi as given 
in (H^l). This yields 

Wiix,x\y,y';gR,fi,D)=TZWi{x,x',y,y';g,D) (3.11) 

where Wi (x, a;', y, y'; g, D) is a regularized expression calculated in D dimensions 
and /z is a subtraction point introduced by the renormalization procedure TZ. 
For the purpose of the present work the only important relation is 

gR = ^,iD-')/^g + o{g^). (3.12) 



While the operation TZ is sufficient to make regular loops well-defined, new 
divergences appear if the contour C has cusps or self-intersections. The renor- 
malization properties of such loops are discussed in [E0| and [Ell. While the 



singularities of regular loops appear at the two-loop level (order 5^ in (3.9)) 
cusps and cross points cause divergences even in leading (non-trivial) order g^. 

Since Wi is indeed plagued by cusps and self-intersections a second renormal- 
ization operation must be carried out to get a renormalized expression W-^ from 
the bare one Wi : According to |^ each cusp is multiplicatively renormalizable 
with a renormalization factor Z depending on the cusp angle. In our case we 
have four cusps with angles 



a 
a' 

P 




= Z{x — x' , —x) , 

= Z(x', X — x') , 

= Ay-y',-y), 

= Ay'-.y-y')- 





= Z{-x,x') 


lyv' 


= A-y,y') 


Ixy 


= A-x,-y) 


Ix'y' 


= A^\y') 


Ix'y 


- A^',-y) 


Ixy' 


= A-^,y') 



(3.13) 
(3.14) 
(3.15) 
(3.16) 



The cross point at the origin introduces a mixing between Wi and 
W2(x,x' ,y,y') :== 

1 tr [C/(0, x) U{x, x') U{x', 0)] ^ tr [C/(0, y) U{y, y') U{y', 0)] ) (.3.17) 

Again the divergences appearing here are functions of the angles 



(3.18) 



The renormalized Wilson loop Wf^ is given by 

W^{x, x', y, y'; gn, ^, Ca, Ca',Cp, C^- , C=j) 

= ^in Z{Ca,gR, fi; D) Z{Ca' , gn, fi; D) Z{Cp, gn, /i; D) ZiCpi , gn, fi; D) 
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X Zii (C^, gn, fi; D) Wi (x, x' , y, y'; qr, fi, D) 
+ Zi2{C^,gB,fi;D) W2{x,x' ,y,y'; gn, fi, D) 



\im^Wi{x,x',y,y';gR,ti,Ca,Ca>,C0,Ci3',C^;D) (3.19) 

where the second renormahzation procedure introduces new subtraction points 
Ca (c.f. and for more details). Of course different renormalization 
procedures are possible and so the Z factors are not unique. We will return to 
this point in Section |^ where we specify a renormalization operation which is 
appropriate for our purposes. 

The observation that Wilson loops with cusps and/or cross points show addi- 
tional divergences has an important consequence for our radial gauge propagator 
as given in ( |3.3| ) : Even the free propagator needs renormalization! This provides 
a natural explanation for the fact that a naive calculation of this object yields 
an ultraviolet divergent result Q . Note that the usual divergences of Wilson 
loops which are removed by TZ, like e.g. vertex divergences, appear at o{g^) and 
thus do not contribute to the free part of the radial gauge propagator, while the 
cusp singularities indeed contribute since they appear at 0(3^) and affect the 



free field case due to the factor 1/g in (3.3) 



Now we are able to answer the question whether the radial gauge is "unphys- 
ical" or "very physical" . It is just its intimate relation to physical, i.e. gauge 
invariant, quantities which makes the gauge propagator — even the free one — 
divergent. One might cast the answer in the following form: The propagator 
diverges because of — and not contrary to — the fact that the radial gauge is 
"very physical". 

Consequently the next questions are: 



— Is there any use for a divergent expression for the free propagator? Espe- 
cially: Can we use it to perform (dimensionally regularized) loop calcula- 
tions? 

— Can one find a renormalization program which yields a finite radial gauge 
propagator? 

In the next Section we will perform a one loop calculation of a Wilson loop using 



the radial gauge propagator (3.7) and compare the dimensionally regularized 
result with a calculation in Feynman gauge. 

In Section ^ we will explicitly demonstrate that the renormalization program 
for link operators carries over and allows to derive a finite result for the radial 
propagator and contrast its properties and use to the regularized version. 
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4 Calculating a Wilson Loop in Radial Gauge 



We choose the path 

(ax , cr e [0, 1] , a; G 

t : z{a) = I w{a -I) , cr e [1, 2] , w{0) = x, w{l) = y (4.1) 
[ (3 - a) y , a e [2, 3] , y e 

It is shown in Fig. |^. The hne w{a — 1) is supposed to be an arbitrary curve 
connecting x and y. 




1 X 
Fig. 3: A Wilson loop with two straight line parts. 



First we will perform the calculation of this Wilson loop in Feynman gauge. 
Using (|3.6|) we get in leading order of the coupling constant 



W{£) 



1 

N 



tr (V exp 



ig f dz^Afj_{z) 



1 + (ig)' 



2N 



3 3 



da dr e{a - r) z^{a) z^t) ^^^^^^""(^(a), z(r)) 




3 3 



,7V^ - 1 1 



2N 2 



l + M'^n^o da dTz^{a)z'^{T)DlT{z{a),z{T)) . (4.2) 







If 



To get rid of the 0-function we have exploited the symmetry property of two- 
point Green functions 



DlT{x,y)^Df^-{y,x). 



(4.3) 
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Decomposing the contour £ according to ( |4.1| ) we find that the Feynman prop- 
agator in (4.2)) connects each segment of £ with itself and with aU the other 
segments. Thus // is given by 



3 3 



(4.4) 



A=l B=l 



where {A, B) denotes the contribution with propagators connecting loop seg- 
ments A and B (c.f. Fig. 0), e.g. 



(1,2) 



1 1 

da /drx^ w^ir) ^'^^"(aa;, w(t)) 







r(D/2- 1) 



da J dr w^{t) [{ax — w{t)) 





211-D/2 



(4.5) 



Next we will evaluate the same Wilson loop in radial gauge. Clearly the first 
and the third part of the path do not contribute if the radial gauge condition 



j_A^^{x) = holds. We insert the free propagator 

{A';,{x)AUy))%^:6'^''D%{x,y) 



(4.6) 



from (3.7) into 



Wi£) = —trlPexp 



ig I da w^(a)A^{w{a)) 



1 + (»5)' 



2N 2 



1 1 

da JdTwf{a)w''{T) Dl^{w{a),w{T)) (4.7) 





and observe that 

Thus the integral in ([4.7j) reduces to 



r(i:>/2- 1) 

47r^/2 



1 1 



da dTw^{a)w'^{T) [{w{a) - w{T)y] 



2ll-£)/2 




1 1 

ds J dt 





-D/2 



(4.8) 
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+ w^(0) «;^(0) [(sw(0) - iu;(0))2]i-^/2 

- w^{l) w^{Q) [{sw{l) - iu;(0))2]i-^/2 

- w^{Q) w^{l) [{sw{0) - 

1 1 

ds /"dTu;^(T)(w^(l)[(su.(l) -w(r))2]i-^/2 







1 1 
dt I da 



-w^(0)[(su;(0) -w(t))2]1-^/2 







w^(0)[(w(a) -iu;(0)) 



211-Z5/2 



(4.9) 



= (2, 2) + (3, 3) + (1, 1) + (3, 1) + (1, 3) + (3, 2) + (1, 2) + (2, 3) + (2, 1) . 



A careful analysis of ( |4.9| ) shows that it exactly coincides with the Feynman 
gauge calculation. This is expressed in the last line where we have denoted 
which parts of the loop are connected by the Feynman gauge propagator to 



reproduce (4.£) term by term. Thus using the radial gauge propagator as given 
in (3.7) yields the same result as the calculation in Feynman gauge. Finally this 
regularized expression has to be renormalized. This can be performed without 
any problems according to ||2^ . Since we are not interested in the Wilson loop 
itself but in the comparison of the results obtained in radial and Feynman gauge, 
we will not calculate the renormalized expression for W{i). 

However a qualitative discussion of the renormalization properties of W{(,) is 
illuminating. By construction W{t) has at least a cusp at the origin. (Other 
cusps are possible at a; or y or along the line parameterized by w, but are not 
important for our considerations.) To give the right behavior of the Wilson 
loop the calculation of W{t) in an arbitrary gauge must reproduce the cusp 
singularity. Usually the parameter integrals in the vicinity of the cusp do the 
job. For gauge choices where the propagator do not vanish in the vicinity of 
the origin this is automatically achieved. Let us assume for a moment that it is 
possible to construct a finite radial gauge propagator obeying the field strength 
formula ( |1.2| ) and therefore have trivial gauge factors along radial lines. Of 
course this is nothing but saying that there are no contributions form parts 
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1 and 3 of the loop, i.e. in the vicinity of the origin. Since the propagator is 
assumed to be finite, there are no singular integrals corresponding to the cusp at 
the origin. Thus a finite radial gauge propagator cannot reproduce the correct 
behavior of the Wilson loop. In turn we conclude that a singular radial gauge 
propagator is mandatory to get the right renormalization properties of Wilson 
loops. 

However, as we will demonstrate in the next Section, the renormalization proce- 
dure for Wilson loops can be used to devise a consistent renormalization program 
for the radial gauges considered here. We will apply it to write down a finite 
version of the free radial propagator. The generalization to higher orders in 
perturbation theory is straightforward. According to our considerations given 
above we shall show that the rcnormalized, thus finite version of the free radial 
propagator is not suitable as an input to perturbative calculations. 



5 The Renormalized Free Propagator 



We define the renormalized radial gauge propagator by 

{A;ix)Aliy)rr''^' lirnJ^^^^^J-^^^-^ lim d;dyW,ix,x',y^y';D) 



y -^y 



(5.1) 

where we have suppressed most of the other variables on which Wi depends (see 



From now on we will concentrate on the calculation of the renormalized free 
propagator {A'^{x)A'^{y))'j^. The details of the renormalization program are 
presented in Appendix Of course the procedure is closely connected to the 
renormalization of cusp singularities of Wilson loops. The result is 

{A';,{x)At{y))% = (5.2) 
= Una (^S'^'f^^-^dy (^^^{n - 7.,) cot 7.,) + (^^(x)At(y))?,d) 

Before discussing some properties of the renormalized free propagator we shall 
show that the counter term 

C;t{x, y) := S'^'ti^-Xdl - 7,,) cot 7.,) (5.3) 



exactly cancels the divergence of the propagator (3/7), i.e. that {A1j^{x)A'l{y)')j^ 
really is finite. To this end we use some tec hnical results derived in Appendix |C|. 
The divergent part of the propagator ( |3.7|) is given by 

1 1 

<<^i^^y) ----^^^^JJ^^^'W jds jdtx-yiisx-tyf]'-^'^ (5.4) 
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6 Summary and Outlook 



Using (|C3) and ( |ai^ we find 



and thus 



1 



r(j/2^i) 



(tt - 7i:y) cot(7r - + finite 

- finite (5.5) 



U^tix,y) + C-^tix,y) = Unite. 



(5.6) 



Note that if one tries to guess a finite expression fike {A'^{x)A^{y))^ one would 
have to introduce a scale /i by hand without interpretation. In our derivation this 
scale appears naturally as the typical renormalization scale of the TZ operation. 

The counter term C^jj{x,y) has some interesting properties. It is symmetric 
with respect to an exchange of all variables and it obeys the gauge condition 



Cll{x,y)=Q = Cll{x,y)y 



(5.7) 



Thus (^°(x)Aj^(2/))5j is finite in the Hmit D — > 4 but still can be interpreted as 
a gluonic two-point function which fulfills the radial gauge condition 



(A«(x)4(y))?j = 0. 



(5.8) 



However the counter term C^'l{x,y) and thus also (A" (x)j4^(y))5j. is ill-defined 
at the origin and hence conflicts with the field strength formula (1.2). Note that 



the regularized propagator in contrast to the renormalized propagator is well 
defi ned and vanishes if one of its arguments approaches zero, as pointed out after 
eq. (3.7). We therefore conclude that we may use the regularized propagator in 
perturbative calculations and can be ensured to preserve relations like the field 
strength formula or eqs. ( 2.11 ) or (3J) throughout the calculation. Although the 
counterterms are not well defined at the reference point itself - a property the 
radial gauge propagator simply inherits from renormalizing the cusp singularity 
of the underlying Wilson line - physical (gauge invariant) quantities are not 
affected, they are rendered finite and unambiguous. 



6 Summary and Outlook 

In this article we have shown how to calculate the radial gauge propagator in a 
Z?-dimensional space using Wilson loops. As discovered in the free propa- 
gator diverges in four-dimensional space. We were able to explain this singular 
behavior by studying the properties of associated Wilson loops. Furthermore we 
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have shown that the free propagator, in spite of being divergent in four dimen- 
sions, can be used for perturbative calculations in a (dimensionally) regularized 
framework and that the result for a gauge invariant quantity agrees with the 
calculation in Feynman gauge. Finally we have presented a renormalization 
procedure for the radial gauge propagator and calculated the explicit form of 
the renormalized free propagator. We have pointed out that any version of the 
radial propagator which is finite in four-dimensional space at least cannot re- 
produce the correct renormalization properties of Wilson loops with cusps at 
the reference point z. 

It is instructive to compare the radial gauge propagators as presented here with 
other approaches: As discussed in Section ^ the radial gauge condition (1.1) 



does not completely fix the gauge degrees of freedom. Thus the field strength 
formula 

1 

Af,{x)^ jdssx^'F.^isx) (6.1) 



is not the only solution of the system of equations^ 

\ F^,{x)^d-A,{x)~d-A^{x). ^''■'> 
One might add a function Q 

Al{x)^dlf{x) (6.3) 



to (3.1) where / is an arbitrary homogeneous function of degree 0. However 
any A^^ [x) added in in order to modify ( |6.1| ) is necessarily singular at the origin. 
Hence regularity at the origin may be used as a uniqueness condition |jll| . If we 
relax this boundary condition other solutions are possible, e.g. 

oc 

A^{x) = - fdssx^F.^isx) (6.4) 



1 



where we must assume that the field strength vanis hes a t infinity. While (6.1) 
is the only solution which is regular at the origin, ( |6.4| ) is regular at infinity. 
Ignoring boundary conditions for the moment one can construct a radial gauge 
propagator by jl4| 

^iG^,{x,y) + G,^iy,x)) (6.5) 

with 

1 CO 

G^,{x,y) := - Jdssx'' Jdtty^ {Fc,^{sx)Fp,{ty)) . (6.6) 

1 



*For simplicity we discuss the QED case here. Aiming at an expression for the free gauge 
propagator this is no restriction of generality. For non-Abelian gauge groups c.f. [ p^ . 
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6 Summary and Outlook 



It turns out that this propagator is finite in four dimensions. However the price 
one has to pay is that boundary conditions are ignored and thus the object 
"hves" in the restricted space \ {0} and not in anymore. In our approach 
we insist on the field strength formula ( |6.lD widely used in operator product 
expansions H and on the regular behavior of vector potentials at the origin 



11 . One might use the propagator (|6.5|) to calculate the ^^-contribution to 



the Wilson loop on the contour (MJJ). It is easy to check that the result differs 



from the one obtained in (4/7), ( |4.9D . Clearly this is due to the fact that (3.5) 
is ill-defined at the origin. 

In the above, all calculations were performed in Euclidean space. In Minkowski 
space Wilson loops show additional divergences if part of the contour coincides 
with the light cone p2[ . Thus we expect the appearance of new singularities also 
for the radial propagator, at least if one or both of its arguments are light-like. 
Further investigation is required to work out the properties of the radial gauge 
propagator in Minkowski space. 

To formulate perturbation theory in a specific gauge the knowledge of the cor- 
rect free propagator is only the first step. In addition one has to check the 
decoupling of Faddeev-Popov ghosts in radial gauge which is suggested by the 
algebraic nature of the gauge condition. However the still continuing discussion 
about temporal and axial gauges might serve as a warning that the decoupling 
of ghosts for algebraic gauge conditions is far from being trivial (c.f. |4| and 
references therein). To prove (or disprove) the decoupling of ghosts in radial 
gauge we expect that our Wilson loop representation of the propagator is of 
great advantage since it yields the possibility to calculate higher loop contri- 
butions in two distinct ways: On the one hand one might use the Wilson loop 
representation to calculate the full radial propagator up to an arbitrary order 
in the coupling constant. The appropriate Wilson loop can be calculated in any 
gauge, e.g. in a covariant gauge. On the other hand the radial propagator might 
be calculated according to Feynman rules. Since the results should coincide this 
might serve as a check for the validity and completeness of a set of radial gauge 
Feynman rules. 
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A Derivation of the Free Radial Propagator 



The free radial propagator derived form the field strength formula shows a 
divergence in Z? = 4, as already indicated in section ||, eq. (3.7). Here we give 
the details of the algebra leading to this conclusion. 



The following relations summarize the steps carried out in the calculation below: 

= |x| , (A.l) 
T^,(x,y) := x'^y^l^g^^d^dy + g^pd^dy-g^^d^d^^-g^f.d^dl) (A.2) 
^ 9t,ud\^\d\y\\x\ \y \ - d^l^ d\y\\y \ - d« yf,d\.^\\x\ + d^d^ x ■ y . 

Introducing x :— x/\x\ and u = s\x\, we have for arbitrary /: 

sxa df'fisx) = x^ d^^fisx) , (A.3) 

1 1 \x\ 

d\x\ jds\x\f{sx) ^ d\x\ jds\x\f{s\x\x) ^ d\x\ Jdu f{ux) = f{\x\x) = f{x) . 



(A.4) 

We get 

{Al{x)Al{y))l^ = 



V{D/2 - 1) ' ^ 
1 1 



jds jdtT^,{sx,ty) [(sx-ty)^] 1-^/2 



1 1 

S-'T.^ix^y) Jds Jdt[{sx~ty)^]'-"^^ 



5"' {g^ud\.,\d\y\\x\ \y\ ~ dlx,d\y\\y\ - dly,d\4x\ + d^dlx- y) 

1 1 

X jds [dt[{sx-tyf]^-^/'^ 







1 1 

- d; Jdsx. [{sx - y)2] 1-^/2 - dyjdty^ [{x - tyf]'-''/' 
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B Renormalization Program for the Free Propagator 



1 1 

+ dpy Jds Jdtx-y[isx~tyf]'-^/^y (A.5) 



1 

The divergent part of the double integral in the last line can be found in Ap- 
pendix At the moment however the exact form of the divergence is not 
important. 



B Renormalization Program for the Free Prop- 
agator 

In Section || we discussed the effect of renormalization on the free radial propa- 
gator. Here derive in detail the appropriate renormalization procedure starting 
form the renormalization properties of Wilson lines. Only a few of the many 
possible renormalization constants will contribute to the final result. 

Since in the relation between the propagator and the appropriate Wilson loop 
( ^.3| ) a factor is involved all quantities especially all the Z's and W's of 



eq. (3.19) have to be calculated up to o{gj^). We have 

m = l + {i9RfSm + o{g^j,) (z = l,2), (B.l) 

Z = l + {tgRfSZ + o{gl,), (B.2) 

Zii = l + iigR)^5Zii+o{g%), (B.3) 

Zi2 = 0+iigR)^SZi2+o{g%) (B.4) 



yielding 



ab 



2N 



Z5-4 



X lim 



D^4 N^-l' 

6Z{C^) + 6Z{Ca') + dZ{Cp) + &Z{C(i,) 
+SZii + SZi2 + SWi 



(B.5) 



Using the fact that up to o((7|.) the two quantities Wi and Wi are essentially 
the same^ we find 



ab 



2N 



lim S ^ , 



lim 5' 



ah 



,D-4 



2N 



lim d'^dy SWi 



hm dldl (1 + {ig)HW^ 



^Only a factor comes in since gji as given in (p.l2[) is dimensionless in contrast to g. 
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lim S' 



ab 



2N 



D^4 - 1 (igy 

= lim(A^(x)4(2/))L. 

U — ^4 



lim d^^Oy Wi 



(B.6) 



To get the SZ^s we must calculate 6Wi which is straightforward using (3.4) and 
( p.ll ). We only need the Feynman propagator (3.6) to get 



SWi 



with 



/|4-_D 



iV2 



ir(D/2- 1) 



-h{x 
-hiv 



,'|4--D 



\y 



2N 47r^/2 

\x-x'r^ + \xf~^ + \y'\^-^ + \y~y'[ 

X — x') + l2{x — x' , —x) + l2{ — X, x') 

y~y') + h{y - y', -y) + h{-y, y') 
-y') + h{x\ -v) + hiy', -x) - hix, -y) 

~x',y- y') + I^{y',~x,y- y') - h{x' -y',x- x',y' - y) 



h 



-y', x-x')+ h{x', -y, X - x') 
1 1 

/i jds Jdte{s-t) 



[(s-<)2]^/2-l 



(B.7) 



(B. 



and 



h{p,q) := ds dt 



P ■ Q 







± 1 

h{'m,p,q) := jds jdt 



[(sp + tq)2]^/2- 



P_<1 

[(m + sp + t(7)2]-D/2-i 



(B.9) 



(B.IO) 



In the following we are interested only in the divergent parts of these integrals. 
The integrals Ii and I2 are calculated in Appendix |^. The results are 



D 



finite 



and 



1 



■ 7 cot 7 + finite 



(B.ll) 



(B.12) 



A- D 

where 7 is the angle between p and q. The integral I3 is finite as long as m ^ 0. 

To specify the rcnormalization factors Z we choose the minimal subtraction 
scheme as described in iQ. In dimensional regularization all the diver- 

gences are given by sums of pole terms. We define every Z factor to be given 
just by the respective sum. The important property of this rcnormalization 
scheme is that the Z factors depend on the angles only and not on the length 
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C Some Important Integrals 



of the loop or of any part of the loop. Using (B.ll and ( B.12D ) the Z factors 
can be read off from (p^) (c.f. 



OZ[Lja) — 


7V2 _ 


1 1 


1 


2N 


47r2 4 


- z? 




- 


1 1 


1 


2N 


47r2 4 


- £1 


SZiCp) = 


N^- 


1 1 


1 


2N 


47r2 4 


- D 


SZ{Cp,) = 


iV2 _ 


1 1 


1 


2N 


47r2 4 


- Z) 


6Zii = 


- 


1 1 


1 


2N 


47r2 4 


- D 


6Zi2 = 


N^- 


1 1 


1 


2N 


47r2 4 


- D 




-{n 


~ Ix'y') 


cot( 



(a cot a — 1) , 
[a' cot a' — 1) , 
(/?cot/3- 1), 
(/?' cot/3' - 1) , 

[(7x'a C0t7j;/y - 1) + (7j;y' COt7a;y/ 

[7j;a;' cot 7x2:' + -fyy' cot 7.yj^' 



(B.13) 
(B.14) 
(B.15) 
(B.16) 

1)] , 
(B.17) 



Now we exploit the fact that only one of the angles, namely ^xyi depends on x 
and y. All the other ones depend only on x or y separately, or on none of them. 
This simplifies ( |B.5| ) drastically: 



{Al{x)Al{y))°j,= ]lm5' 



ab 



2N 



D-4 



D^4 N^-l' 

/I 1 



lim d^dy 5Zi2 + 5W, 



(B.19) 



= hm. ( 5'^'^i''-^^l^l ( —^^{ir - l.y) cot 7., + {Al{x)Al{y))% 



where we have used (B.6) to get the last expression. 



C Some Important Integrals 



The integrals Ii and I2 played an important part in the renormalization pro- 
cedure of appendix ^ and determine the divergences of the naive free radial 
propagator introduced in section ||. They are discussed in detail below. 



To calculate 



1 1 

Ji := jds Jdte{s-t) 





we introduce the substitution 



[(5-i)2]W2-l 



(C.l) 



g = s~t , h = s + t (C.2) 
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to get 



1 2-3 1 

1 Jdg jdhg^-^^ jdg{l-g)g'-^ 



2 

3 

r(2)r(3-L') 1 



(C.3) 



For the calculation of 



1 1 





we have to distinguish the two cases p aq where the only divergence that 
appears is for s ^ t = and p — aq with an additional divergence at s = ta. 
Here we will only need the former. 

As a first step it is useful to separate off the divergence at the origin by the 
substitution 

X = s + t , x^s/X. (C.5) 



This yields 



1/2 l/{l-x} 1 1/x 



l2{p,q) = \ Idx I dX+ Idx IdXlX^-'^ ^'^ 



1/2 



[{xp+il-x)qy]D/ 



2-1 



4-D [{xp+{l-x)q)^]Df^- 



1 



1 

D— 4 

, X P_q . 

A-D [{xp + (1 - x)g)2]0/2-i ■ ^^■^> 

1/2 

As long as p ^ —q holds there are no divergences in the a:-integration since 

u[x) -.^ xp + [I ~ x)q (C.7) 
never vanishes. We introduce the angle between p and q 

cos 7 (C.8) 

\p\ \q\ 

and the substitution [|8) 

e-'^ = + . (C.9) 

x\p\ + {1 - x)\q\e-'i ^ ' 
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C Some Important Integrals 




Fig. 4: The geometry of the variables appearing in the calculation of I2 (C.4) 



Note that -0 is nothing but the angle between p and u{x). To perform this 
substitution in ( p.6D we need 

(C.IO) 
(C.ll) 



X — \q\ sin{-f — ip) / N {tp) , 1 — x = \p\sm^/N{^/j) , 



uix)Y = p^q^ sin^ -f/[N{i}})Y and -- = - 



dx \p\ \q\ sin 7 



with 



N[il}) := \p\ saup + \q\ sin(7 — ip) . 
In addition it is useful to introduce 

i)' := i}j{x = 1/2) 

which is the angle between p and p + q (cf. Fig. %. 

Using all that we end up with 



(C.12) 
(C.13) 



h{p,q) ^ dip 



-\p\ \q\ sin7 / \p\ sinip 



7V2 



N 



D-4 



9 9 • 2 
p'^q'^ sm 7 



D/2-1 



p-q 



d^/j- 



\p\ \q\smj f\q\ sm{j -tp) 



iV2 



N 



D-4 



4-D 
^2 X D/2-1 



p^q^ siv? 7 



P-q 

A- D 



cos 7 sin'^ ^ 7 
4-D 



cos 7 sin'^ ^ 7 







\q\'''^ jdiP sin^-^V' + \p\''-^ jd^P sin^-^(7 - V') 



|gr^yd^sin^-4i^ + |p|4-^ y d?Asin^-4^ 
7 7 



4-i:) 



7 cot 7 + finite . 



(C.14) 
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